Critical dynamics and effective exponents of magnets with extended impurities 
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We investigate the asymptotic and effective static and dynamic critical behavior of (d — 3)- 
dimensional magnets with quenched extended defects, correlated in Ed dimensions (which can be 
considered as the dimensionality of the defects) and randomly distributed in the remaining d — Ed 
dimensions. The field-theoretical renormalization group perturbative expansions being evaluated 
naively do not allow for the reliable numerical data. We apply the Chisholm-Borel resummation 
technique to restore convergence of the two-loop expansions and report the numerical values of the 
asymptotic critical exponents for the model A dynamics. We discuss different scenarios for static 
and dynamic effective critical behavior and give values for corresponding non-universal exponents. 

PACS numbers: 61.43.-j, 64.60.Ak, 75.10.Hk 



I. INTRODUCTION 



Critical properties of structurally disordered magnets 
remain a problem of great interest in condensed mat- 
ter physics, since real magnetic crystals are usually non 
ideal. A simple and natural case of disorder is imple- 
mented via the point-like uncorrelated quenched non- 
magnetic impurities and is experimentally realized as 
substitutional disorder in uniaxial Q as well as in Heisen- 
berg 0, Q magnets. Examples are given by substitute 
alloys Mn x Zni_ x F2, Fe x Zni_ x F2 for the uniaxial (Ising) 
magnets Q, and by amorphous magnets Fe9o+ x Zrio- x , 
Fego- y M y Zrio ( M = Co, Mn, Ni) Ull, transition-metal 
based magnetic glasses 0, as well as disordered crys- 
talline materials Feioo- x Pt x 9j, Fe7oNi3o |8| and Eu- 
chalcogenide solid solutions for the Heiscnberg mag- 
nets. The question of a great interest arising here is: does 
the disorder change critical properties of the systems? 
The answer is given by the famous Harris criterion 
It states, that disorder changes the critical exponents, if 
the critical exponent a p of the pure system is positive: 
a p = 2 — dvp > 0. Here, v v is the critical exponent 



governing the divergence of correlation length of the cor- 
responding pure system and d is the space dimension. 

On the base of this inequality, one can estimate the 
marginal value m c for the spin m-vector model of mag- 
netic systems, such that for to > m c the critical expo- 
nents remain unchanged by point-like defects whereas for 
?7i < m c they cross over to new values. The present esti- 
mates for m c in three dimensions definitely imply < 2: 
m c = 1.942 ± 0.026 [13 and m c = 1.912 ± 0.004 [3 and 
thus only the pure Ising model (m = 1) is affected by 
weak point-like uncorrelated disorder at criticality. 
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But in real magnets one encounters non-idealities of 
structure, which can not be modelled by simple point-like 
uncorrelated impurities. Indeed, magnetic crystals often 
contain defects of a more complex structure: linear dis- 
locations, disclinations, complexes of non-magnetic im- 
purities, embedded in the matrix of the original crystal 
|14| . Theoretical studies of critical behavior of magnets 
containing such "extended" (macrosc opic) d efects have 
attracted considerable interest [TH ITfl UtI \1§L Il9l I20L I2H 
I22I E^L I24L Efil . EfiT ] , however, a systematic experimental 
analysis still remains to be performed. 

One of the possible ways to treat defects that extend 
through the system being randomly distributed in space 
but oriented in the same direction is to consider them as 
quenched e^-dimensional nonmagnetic impurities of par- 
allel orientation. This was proposed in the work of Doro- 
govtsev . The case Ed = is associated with point-like 
defects, and extended parallel linear (planar) defects are 
described by 6d — 1(2). To give an interpretation to the 
non-integer values of Ed, one may consider patterns of 
extended defects like aggregation clusters, and treat Ed 
as the fractal dimension of these clusters [ljj- However, 
relation of analytically continued non-integer Euclidean 
dimension to the fractal dimension is not straightforward 

M 

It was shown |16|. that presence of extended impuri- 
ties leads to a generalized Harris criterion, namely, in 
this case disorder alters the critical behavior of the pure 
system, if Ed > d — 2/v p . If one considers the point- 
like disorder with Ed — 0, the generalized Harris criterion 
turns to ordinary one cited above. Again, this inequality 
defines for each value of Ed the critical value m c , below 
which the extended disorder alters the universality class. 
But now the disorder induced critical behavior of d = 3 
systems holds not only for the Ising magnets (m = 1) but 
for m > 1 as well, as shown in Fig. ^ (region denoted 
as "Diluted" in the figure). Therefore, the class of mag- 
nets where the new critical behavior may be found is not 
restricted to the uniaxial ones (m =1), but include easy- 
plane (m = 2) and Heisenberg (to = 3) systems. Thus 
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predicted phenomena of a new universal behavior may 
be experimentally checked for a wider class of magnets. 




FIG. 1: The phase diagram of d = 3-dimensional m- 
vector magnets in the presence of Ed- dimensional impurities. 
Marginal value m c (ed) separates regions of different univer- 
sality classes. Disorder is relevant in the £d-m plane below the 
m c curve. Filled squares: results of calculations, based on six- 
loop renormalization group data 28j for the correlation length 
critical exponent, dotted line: estimate based on resummed 
e, ed-expansion for m c , Eq. 13UII . Only the physical region 
£d > is shown. See Section ITVI for a whole description. 

Another interesting feature of systems with parallel ex- 
tended defects is that due to the spatial anisotropy they 
are described by two correlation lengths, one perpendic- 
ular, and one parallel, to the extended impu- 
rities direction ^5|. As the critical temperature T c is 
approached, their divergences are characterized by corre- 
sponding critical exponents v\ \ : 



t 



in 



1*1" 



(1) 



where t is the reduced distance to the critical tempera- 
ture t = (T — T c )/T c , Anisotropic scaling holds also for 
the spin-spin pair correlation function at T c and is gov- 
erned by exponents r]± and n \ . Whereas the magnetic 
susceptibility is isotropic, as far as all order parameter 
components interact with defects in a similar way. How- 
ever, the dynamic critical behavior is modified, two times 
of relaxation in directions perpendicular and parallel to 
the extended impurities t± and ti behave correspond- 
ingly as: 



1~L 



Si 



(2) 



with dynamical exponents z± and z\ \ . 

Therefore, magnets with quenched extended defects of 
parallel orientation constitute a large class of systems 
with a bulk of unusual phenomena worth to be ana- 
lyzed. Their asymptotic critical behavior was a subject of 



the field-theoretical renormalization group (RG) analyses 
[3 H [H IH E3 . In particular, a double expansion in 
both e = 4 — d, Ed was suggested and RG functions were 
calculated 0] to order e, Ed] qualitatively, the crossover 
to a new universality class in the presence of extended 
defects was supported. These calculations were extended 
to the second order in Refs. (1(3 . [l8| . However, these 
divergent RG expansion did not give a reliable numerical 
estimate for the critical exponents - the goal highly de- 
sirable both for experimental and simulational purposes. 
Numerical estimates for the exponents describing static 
critical behavior of d = 3 magnets with parallel extended 
impurities were obtained only recently |2fJ, |21j by apply- 
ing special resummation technique to the two-loop RG 
functions of Refs. fl6l [l8|. Currently, no estimates of 
similar accuracy exist for the dynamic exponents describ- 
ing critical slowing down in the magnets under consider- 
ation. The theoretical investigation of critical dynamics 
of system with parallel extended defects was performed 
in the one-loop approximation |17| for the simple models 
with non-conserved order parameter (model A) and con- 
served order parameter (model B). Then the RG study 
of model A critical dynamics was extended to the two- 
loop order 01 • But again, the convergence properties of 
the series obtained did not allow numerical evaluations. 
Recently, the short-time critical dynamics of the model 
A was considered in Ref. 1221. 



The goal of this paper is to apply the state-of-the-art 
analysis of the divergent RG perturbation theory series 
to get numerical estimates of the exponents describing 
model A dynamics of d — 3 magnets with extended im- 
purities for a wide region of impurity dimension Ed and for 
different values of to. Complementing existing estimates 
of the static exponents |2(J, |21| , our results will give com- 
prehensive values for the critical exponents of systems 
with extended impurities and, in this way, will facilitate 
their experimental analysis. Moreover, since both in ex- 
periments and in Monte Carlo (MC) simulations often 
the asymptotic region is not reached and effective expo- 
nents are observed [2{j, we will calculate these as well. 
By these estimates we will predict the possible scenarios 
of approaching the critical point. These are important 
to perform corresponding experimental checks, as it was 
shown recently by theoretical [3(j and experimental [3l| 
studies of the d = 3 disordered Heisenberg magnets with 
point-like defects. 



The setup of the paper is as follows. In the next sec- 
tion we present the model; in section ITTT1 the renormaliza- 
tion procedure is discussed. In section llVl we apply the 
resummation techniques to analyze the renormalization 
group functions in two-loop approximation and present 
the quantitative estimates for the asymptotic dynamical 
critical exponents. SectionlVl gives description of possible 
scenarios for the effective critical behavior and section IVT1 
concludes our study. 
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II. THE MODEL 



The starting point is the effective Hamiltonian of the 
modei of an m-vector magnet with e^-dimensional de- 
fects, extending throughout the system along the coor- 
dinate directions symbolized as x\\ and randomly dis- 
tributed in perpendicular directions x± |!6lll8j| : 



H = J d d x[-(^ + V(x))^(x) + (Vj(x)f 
+ ao (V||^)) 2 + ^(0 2 (x)) 2 ]. 



(3) 



correlation: 

(rf{x,t)r?{x',t')) = 2\ 5(x-x , )6(t-t , )6 lj . (6) 

The brackets (...) mean an average with respect to the 
thermal noise. 

Within the field theory approach it is convenient to 
use the Bausch-Janssen- Wagner formulation |3i| which 
is given by the Lagrangian: 



d d x dt ^2 ' 



dcj) 1 8H 



(7) 



Here, (j) is an m-component vector field: <fi = {4> l ■ ■ ■ </>™ 1 }, 
fio and uq are the bare mass and the coupling, ao is the 
bare anisotropy constant, V|| and V± denote differen- 
tiation in the coordinates x\\ and xj_ and the impurity 
potential V(x) is introduced. The probability distribu- 
tion of defects has zero mean and variance given by: 



((V(x)V(y))) 



v 6 d £d (xj_ -yx). 



(4) 



Here, ((...)) denotes the averaging over the distribution 
of defects, (-«o) is a positive coupling constant, which 
is proportional to the concentration of impurities. The 
constant ao parameterizes the anisotropy, arising in the 
system due to the presence of extended defects. 

We consider critical dynamics of the model © for the 
case of a non-conserved order parameter. For this case 
the d yna mics can be expressed in the Langevin equation 
form 321: 



dcj) l (x,t) &H 
— — Ao- 



Ot 



d(f) l (x,t) 



+ rf(x,t), i = l...m, (5) 



where A is the Onsager kinetic coefficient and r] l (x,t) is 
the Gaussian random-noise source with zero mean and 



Here, <j> 1 are components of an auxiliary response field 
introduced to average over the thermal noise. Then cor- 
relation and response functions are computed with the 
help of a weight function W ~ q-^M . 

It is known that studying dynamical properties of dis- 
ordered systems averaging over random impurities can 
be applied directly to dynamical weight function W. As 
it is established in Ref. j3^|, the configurational averag- 
ing can be performed avoiding replica trick |35|. How- 
ever it leads to the same perturbative expansions for the 
RG functions as those obtained with replica formalism. 
Therefore both approaches giving equivalent results are 
possible for dynamics. However the above alternative 
does not exist when the static critical behavior is ana- 
lyzed (see e.g. 0,0). As far as both static and dynamic 
criticality is addressed in this paper, we use the replica 
trick, representing the logarithm of a weight function in 
the following form: 



In W = lim 



((W n - 1» 



(8) 



Finally it leads to study of properties of replicated La- 
grangian £[4>,4>] [3(|: 



C[$,4>] = W J d d xdt Y,fta 



^ + A (/i 2 - Vi - ooVf )C - A O 0L + E 



E A oy / ^d rf tf«ftdt , *(a: ± -vj.)^(i > *)^(a :> *)^(i/ J t / )^(l/ J < / )k 



(9) 



here, the summation over Greek indices spans from 1 to 
n denoting the different replicas and the Latin indices go 
from 1 to m denoting the components of the order param- 
eter. To study behavior of this model J5| in the vicinity 
of the critical point we apply the minimal subtraction 
scheme within field-theoretical RG. The description of 
this approach is given in the next section. 



III. RG STUDY 



The description of the long-distance properties of the 
model © near the second order phase transition point 
is performed using the field-theoretical RG method [3^| . 
Let us present the renormalization algorithm, developed 
for Lagrangian field theory. It is well known, that such 
a theory encounters with ultraviolet divergences, the re- 
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moval of which is achieved within an appropriate renor- 
malization procedure by a controlled rearangement of the 
perturbation theory series. The change of bare couplings 
Mo, Vq and of the anisotropy constant ao under renormal- 
ization is described by the RG functions: 



p u {u,v) 

Ca(u,v) 



du 

d In k 
d In a 



|o, (3 v (u,v) = 



dv 
dint 



o- 



(10) 
(11) 



Here, u, v, a are the renormalized couplings and 
anisotropy constant, respectively, k is the rescaling pa- 
rameter, the notation |o indicates differentiation at fixed 
bare parameters. The bare fields <p, (f>, the mass [iq and 
the Onsager kinetic coefficient Ao are related to the renor- 
malized ones ip, (p fx and A by: 



2 rr 2 

MO = Z M 2 M ) 



A 1 = Z\\ 



where Z-factors are dimensionless functions of renormal- 
ized parameters m,u,v,a. Their flows are defined by 
corresponding RG functions: 



C(u,v) 



d\nZ v 
d In k 

dlnn 
dbxZ~i 

d In k 



70, 



d\nn 

i are connected by the relation: 



(12) 
(13) 

(14) 
(15) 

C = 



Functions £ and 

(70 " 7*)/2- 

The fixed points (FPs) u*,v* of the RG transformation 
are defined as common zeroes of the /3-functions: 



(3 u {u*,v*) = 0, /?»*) = 0. 



(16) 



A FP is stable, if the eigenvalues of the stability matrix, 
defined as: By = d(3 Ui (u* , v*)/d(3 U] (u*, v*), m = {u,v} 
have positive real parts. In the case when the stable FP 
is physically accessible, i.e. can be reached starting from 
the initial values of the renormalized couplings (in our 
case, u > 0, v < 0), corresponds to the critical point of 
the system. The RG functions (|12[) - (|15fl taken at this 
point give the critical exponents of magnetic susceptibil- 
ity, correlation length and relaxation time: 



1 



7^2 (u*,v*) 



2 - 70 (u*, v*) ' 
2 - 702(u*,u*) - 70(u*,i;*), 
2 + ((u*,v*). 



(17) 
(18) 



Note, that the following relations between the expo- 
nents describing the parallel and perpendicular correla- 



tion length and relaxation times hold |15| : 
*|| = W(l- y)- 



(19) 



To obtain the quantitative characteristics of the dy- 
namical critical behavior of magnetic systems with ex- 
tended impurities, we turn our attention to the RG func- 
tions derived in Ref. |18| in two-loop approximation: 



Pu/u = -£ + 



12 



(m + 8) n (3m + 14) 2 
■u + 2v- - — -u 1 



6 



12 



g(llm + 58)+(m-4) ^ 
6 6(e + Ed) 



144 



328 + 32- 



£d 



Pv/v = 

70 = 
70 2 



-e - s d + -v 



£ + £d 
m + 2 



m + 2 



11 - 



u v 

3 6 



36 36 
m + 2 1 



£ + £d 
-vu + 



5 m + 2 

? 

12 3 



u- 



(a 



6 

m + 2 

24 
1 5 



r- 



72 

m + 2 



G 



u 2 - 2Av 2 



-vu 



6 - 



£d 



£ + £d 



2 {m + 2) 

— V V vu, 

3 36 36 



1 (m + 2)(61n 

v H 

3 72 

(m + 2) 5 2 



1) 



36 



vu H v . 

36 



(20) 

(21) 
(22) 

(23) 
(24) 

(25) 



Here, e = 4 — d and the replica limit n = has been 
assumed. Putting in Eqs. (20)1 -(SO), (ESJ £d = one 
regains the RG functions of the m-vector magnet with 
point-like defects pll39ll40l|. 

The description of the critical behavior of mag nets 
with extended impurities carried out in Refs. |l5l Il6l 
ITsl IT^ | was based on the double expansion in two param- 
eters £, £d, considering both to be of the same order of 
magnitude. Being quite clear technically, such a state- 
ment causes however certain cautions. Indeed, taking 
the e as a small parameter is due to the fact that it de- 
viates from the upper critical dimension, where the (/> 4 
theory becomes asymptotically free |37). whereas £d is 
the dimension of defect itself and obviously has a differ- 
ent physical origin. Therefore, it is desirable to search for 
an alternative way to analyze the RG functions (|20l) - (|25|l . 
Fortunately, such an alternative exists and it is exploited 
in the field-theoretical approach to critical phenomena 
[4l| . Namely, it consists in the analysis of minimal sub- 
traction RG functions directly at dimension of space of 
interest, fixing e and solving the FP equations numeri- 
cally [IJ (the so-called fixed dimension scheme). It was 
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proposed |2(J, |2l| to extend the approach of direct eval- 
uation to the RG functions of the present model, i.e. to 
treat them directly at d = 3 (e = 1) for different fixed 
values of the (non- integer) defect dimensionality Ed- In 
our calculations (see Section IIV|) , we will make use of 
both possibilities, exploiting e, Ed expansion as well as 
working within the fixed dimension scheme. 

In the field-theoretical RG approach the series expan- 
sions of the RG functions in powers of the coupling ap- 
pear to be divergent; moreover, they are characterized by 
a factorial growth of the coefficients implying a zero ra- 
dius of convergence [32} • To take into account the higher 
order contributions, the application of special tools of 
resummation is required |42j. In our previous paper 
|2l|. analyzing static critical behavior we applied the 
Chisholm-Borel resummation technique. Here, the Borcl 
image of the initial function is extrapolated by a rational 
Chisholm |43| approximant [K/L}(x, y). One constructs 
this ratio of two polynomials of order K and L such that 
its truncated Taylor expansion is equal to that of the 
Borel image of the initial function. The resummed func- 
tion is then calculated by an inverse Borel transform of 
this approximant. The details can be found in |44| . 



IV. THE RESULTS 



critical behavior is expected in the region of Ed — m plane 
denoted as "Diluted" . In the region denoted "Pure" the 
asymptotic critical behavior does not indicate presence of 
defects, however the effective critical behavior does. Be- 
low, we will perform an analysis of the critical dynamics 
for both regions of the phase diagram. 

A. Fixed d = 3 scheme 

We start from the RG /^-functions l(2Tjj) , (|2*T)) : fixing the 
value e = 1 (i.e. d = 3) and treating Ed as a varying pa- 
rameter, we look for the common zeros of the resummed 
functions (3 U and (3 V . The numerical values of the sta- 
ble FP coordinates obtained for the three-dimensional 
m-component magnets with m = 1,2,3,4 can be found 
in our previous paper |2l|. 

To calculate the values of the critical exponents z± we 
substitute Eq. (|25|) for £ in Eq. (|18|l . apply the resum- 
mation procedure for the resulting series, and, finally, 
estimate them at the stable FP. Note, that because the 
function ^ is not symmetric in variables u and v (namely, 
it does not contain term linear in u) the Chisholm approx- 
imant chosen for its analytic continuation differs from 
those chosen for the (3- functions. The value of critical 
exponent z\\ is obtained using the relation 1191) . 



The critical behavior of the model is influenced 
by the presence of the following FPs: the Gaussian FP 
G (u* = v* =0), the 0(m)-symmetric FP of a pure 
magnet P (u* ^ 0,v* — 0) and the random FP R 
(u* 7^ 0, v* 7^ 0) that governs disorder-induced critical 
behavior. A polymer FP with u* — 0, v* > is not 
reachable from the initial values of couplings and there- 
fore is out of interest for our study. Depending on the 
values of global parameters d,£d,m one of the above FPs 
is stable and in the asymptotic limit governs the crit- 
icality. However, as we will show in the section an 
approach to the asymptotic limit and hence the effective 
critical behavior is influenced by all the FPs. 

For the d — 3 magnets considered here, it is the 
crossover between FPs P and R that corresponds to 
the change of the universal properties of an m-vector 
magnet upon dilution by e^-dimensional defects. Ac- 
cording to the generalized Harris criterion, this crossover 
occurs at certain marginal value m c (£d). At m > m c , 
the FP P is stable, indicating no change in critical be- 
havior, whereas for m < m c one finds the stability of 
the FP R, displaying the fact of relevance of disorder. 
From the generalized Harris criterion one can obtain 
the marginal value £™ arg as function of m. Using the 
six-loop results j2^| for the correlation length critical 
exponent v p (m) of the pure m-vector magnet, one ob- 
tains: e™ aiB (m = 1) = -0.173; e™ aie (m = 2) = 0.016; 
e™ arg ( TO = 3) = 0.172; e^{m = 4) = 0.300. These es- 
timates are shown in Fig. ^by filled squares connected by 
a solid line. The figure may serve as a phase diagram of 
an m-vector magnet with extended impurities: the new 



TABLE I: The values of dynamical critical exponents of 
three-dimensional m-component magnets at different fixed 
values of extended defect dimension Ed- 

m = 1 m = 2 m = 3 m — 4 



Sd z± 

2.172 

0.1 2.248 

0.2 2.302 

0.3 2.340 

0.4 2.367 

0.5 2.386 

0.6, 2.399 

0.7 2.408 

0.8 2.413 

0.9 2.416 

1.0 2.418 

1.1 2.418, 



2.094 
2.111 
2.127 
2.142 
2.156 
2.169 
2.181 
2.194 
2.206 
2.217 
2.228 



z± 
2.065 
2.079 
2.135 
2.183 
2.223 
2.256 
2.284 
2.307 
2.326 
2.342 
2.356 
2.368 



2.068 
2.073 
2.074 
2.081 
2.087 
2.092 
2.096 
2.100 
2.104 
2.107 
2.109 



2.062 
2.062 
2.062 
2.095 
2.133 
2.166 
2.195 
2.219 
2.241 
2.260 
2.276 
2.290 



z \\ 

2.062 
2.062 
2.063 
2.063 
2.064 
2.067 
2.069 
2.070 
2.072 
2.073 
2.073 



z± 
2.058 
2.058 
2.058 
2.058 
2.083 
2.111 
2.138 
2.162 
2.184 
2.203 
2.219 
2.233 



2.058 
2.058 
2.058 
2.061 
2.062 
2.065 
2.067 
2.069 
2.070 
2.071 
2.071 



In Tabled we give the obtained results for critical ex- 
ponents zu, z± of three-dimensional m-component mag- 
nets. The case Ed = corresponds to point-like defects. 
As it was already noted in the Introduction, for this case 
only the Ising model (m = 1) is influenced by disorder, 
and thus for systems with m > 1 the values of the criti- 
cal exponents are not altered by disorder. Here, one can 
compare our result with the dynamic exponent of the 
pure and diluted Ising models, as shown in Tables ITTllIIII 
When Ed increases, for m = 2, 3, 4 the critical exponents 
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remain constant and equal to the corresponding expo- 
nents of the pure model, until m becomes m c for given 
£d and for m > m c the values of exponents start to in- 
crease because they take their new values belonging to 
the disordered universality class. 

TABLE II: The value of dynamical critical exponent z for 
pure three-dimensional Ising model. MC: Monte Carlo sim- 
ulations; exp: experimental estimates for FeF2. Theoretical 
estimates, ~ e 2 (e 3 ): direct substitution e = 1 into ~ e 2 (e 3 )- 
expansions; ~ e 3 , res: resummation of e 3 expansion; 4 loops, 
res: resummation of four-loop massive 3d RG functions. 





Method 


Result 


Ref. [4£] 


MC 


2.04 ± 0.03 


Ref. [4JJ 


MC 


2.06 ±0.02 


Ref. [48] 


MC 


2.03 ±0.01 


Ref. [49] 


MC 


2.032 ± 0.004 


Ref. [£0J 


MC 


2.055 ± 0.003 


Ref. [51] 


exp 


2.1 ±0.1 


Ref. [18] 


~e 2 


2.014 


Eq. 


~e 3 


2.024 


Eq. 


~ e 3 , res 


2.012 


Ref. [52] 


4 loops, res 


2.017 



TABLE III: The value of dynamical critical exponent z for 
three-dimensional Ising model with point-like uncorrelated 
defects. MC: Monte Carlo simulations; exp: experimental 
estimates for Feo.9Zno.1F2. Theoretical estimates, ~ e 1//2 : di- 
rect substitution e = 1 into e 1//2 -expansions; 2 loops, res: re- 
summation of two-loop massive 3d RG functions; ~ e 1/2 (e 3/2 ) 
res: substitution of the two (three)-loop resummed FP coor- 
dinates into the two- loop expansions for z. 





Method 


Result 


Ref. [53] 


MC 


2.13 ±0.05 


Ref. [54] 


MC 


2.4 ±0.1 


Ref. [£5J 


MC 


2.62 ±0.07 


Ref. [52] 


exp 


2.18 ±0.10 


Ref. [57] 


,-1/2 


2.336 


Ref. [5J 


2 loops, res 


2.237 


Ref. [38] 


~ e 1//2 , res 


2.023 


Ref. [58] 


~ e 3/2 , res 


2.180 



An interesting feature of the data shown in Table [I] is 
that for each fixed Ed the relation zu < z± holds. One 
can give the following physical interpretation to this fact. 
As it has been noted for isotropic systems j32, the dy- 
namical critical exponent z is proportional to the pair 



correlation function critical exponents 77: 

z = 2 + cr). (26) 
For the systems we consider here the critical exponents 
T7_i_ and t]\ \ , that characterize the behavior of the pair spin- 
spin correlation function in the directions, perpendicular 
and parallel to the extended defects, are distinguished 
(lif . As far as the extended defects cut interacting bonds 
of spins perpendicular to the extended-defect direction, 
in the parallel direction the fluctuations are stronger and 
therefore 7711 < i]± and thus, by Eq. (|26[1 with 3511, z± 
we obtain the confirmation to our results. However, in 
the anisotropic case the coefficient c in Eq. I|26[) is also 
direction dependent. Thus the above argumentation is 
rather of qualitative nature to give a hint to physical 
interpretation of the observed relation zu < z±_ . 

In our previous papers [13, El we have also touched 
the question of existence of an upper marginal value for 
the defect dimensionality Ed- In our analysis, we observe 
the disappearance of a stable reachable FP for Ed slightly 
above 1. The following physical interpretation was pro- 
posed: extended defects of large dimension (e.g. paral- 
lel planar defects with Ed = 2), may divide the system 
into non-interacting regions and thus inhibit ferromag- 
netic order. Recently in Refs. 0] the Ising magnet with 
spin interaction bonds, correlated in 2 dimensions and 
randomly distributed in the perpendicular direction, was 
studied. Although such a system differs from those con- 
sidered here, both models possesses a number of common 
features. In particular, the smearing of a phase transition 
due to the presence of planar defects was predicted and 
explained by the existence - although rare - of infinite 
spatial regions, which are free of defects and therefore 
may be locally in the ordered phase. 

To confirm the obtained results, we have also tried to 
apply the Pade-Borel resummation to estimate the val- 
ues for z±, Zu. In this procedure, one rewrites the two- 
variable series (J20J), i|21|l in terms of a resolvent series 
|59| of one variable and then applies a Pade approximant 
for its analytic continuation. We do not present the re- 
sults obtained, they give rather close numbers to those 
in Table ITTT1 and reproduce the same behavior. 



B. e, Ed-expansion 

Another possibility to obtain estimates for dynamic 
critical exponents of the m-vector magnet with extended 
defects is to resum the double s, Sd ex pan sions obtained 
in the two- loop approximation in Ref. |18| . For m < m c , 
the random FP R is stable and the expression for z±_ and 
anisotropy function £ a for m/1 read 0] : 
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z± = 2 - A ™ + \ , e + j. m + V. g + {-4(m + 2) [5m 2 + 42m + 112 - 192 (m - 1) ln(4/3)]e 2 
4(m— 1) 8(m — 1) 

- 4(m + 2) [27m 2 - 264m - 240 + 576(m - 1) ln(4/3)]ee 

+ [59m 3 - 528m 2 - 2928m - 896 + 1728(m + 2)(m - 1) ln(4/3)]e 2 }(1024(m - l) 3 ) -1 ; (27) 
(m + 2) m + 8 _ 
Ca = 4(m- 1) £+ 8(m- 1) £ 

+ [-4(m + 2) (5m 2 + 10m + 144)e 2 - 4(m + 2) (27m 2 - 168m - 336)ee 

+ (59m 3 - 240m 2 - 2640m- 1472)e 2 ] [1024(m - l) 3 ] -1 , e = e + e d . (28) 



Whereas for m > m c , the pure FP P is stable. There- 
fore the critical behavior is isotropic and the dynamic 
exponent z reads: 

z = 2 + 0.363 ( TO + 2 )f . (29 ) 
(m + 8) 2 

As it is known 60], due to degeneracy of the (3- 
functions of the weakly diluted Ising model with point- 
like defects, the usual e-expansion for critical exponents 
turns into the -y^e-expansion. The same holds for the ex- 
tended defects: indeed, the case m = 1 is to be analyzed 
separately as far as expressions (|27|l . (|28|l contain poles 
at m = 1, and one arrives [Ts[ at the -^/e-expansion for 



the critical exponents. In particular, the corresponding 
expressions for z can be given only to the lowest non- 
trivial order. Moreover, the i/e-expansion does not allow 
for a reliable numerical estimate |4f^ . Therefore, for the 
disordered Ising model, the fixed dimension scheme con- 
sidered in the subsection IIV Al remains the only way to 
get numerical estimates. 

To get the numerical estimate for m c in the frames 
of e, £d-expansion we substitute into the modified Harris 
criterion the five-loop e - expansion for the correlation 
length critical exponent v v for pure m-component model, 
given in [6lJ and obtain the following expansion: 



fee 2 e 3 e 4 
m c = -4 + 8- - 2.50000 1.500000s - 2.448919— + 3.014557e 2 + 4.141561ee - 1.682130 

\ s e e e 

- 14.12940 e 2 i - 0.5736055 e e 2 + 7.657623 s 3 + 55.57104 s 3 e - 16.25104 e 2 e 2 - 37.62878 e e 3 + 

+ 22.53257 e 4 - 3.345417 7) ( 2 " ^) > e = e + e d . (30) 



Putting here e = e (i.e Sd = 0) one recovers the e- 
expansion for m c of the model with point-like uncorre- 
cted defects 

To estimate m c numerically for different fixed values 
of <i,£rf, one should apply a resummation. As it is known 
0,|62j already simple Pade-analysis gives convergent re- 
sults for marginal dimensions. In the Pade analysis [&|, 
it is known that the best convergence is achieved along 
the main diagonal of a Pade table, therefore we make use 
of the diagonal [2/2] Pade-approximant. In the three- 
dimensional case (e = 1) it gives: m c (ed — 0) = 1.92 
(which is in a good agreement with the known six-loop 
results for the point-like disorder m c = 1.942 ± 0.026 
[12 and m c = 1.912 ± 0.004 p]J), m c (e d = 0.1) = 2.48, 
m c (e d = 0.2) = 3.10, m c (e d = 0.3) = 3.77, m c (e d = 
0.4) = 4.54. The results obtained are plotted by the 
dotted line in Fig. ^ 



To estimate dynamical critical exponents in the FP R 
we applied the Chisholm-Borel resummation technique 
to the functions Ij27(l . Ij28(l treating them as series in two 
variables e,Ed and using the Chisholm approximant. In 
the corresponding two-loop approximation the dynamical 
critical exponents in the pure FP P is given by Eq. Ij29(l 
and the expression is too short to be resummed. There- 
fore it has been estimated at e = 1 and different m by 
a direct substitution. The results are presented in Ta- 
ble [W] They qualitatively confirm our conclusions made 
on the basis of fixed dimension technique in subsection 
IIV Al (see Table llllll . The discrepancy between the data 
of the two techniques may serve also as an estimate for 
accuracy of the perturbation schemes applied. 

Let us note, that currently the three-loop series for 
the exponent z of the pure m-vector magnet model A 
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dynamics is available. The expression found in 6J| reads: 
z = 2 + 0.726(1 -£ -0.1885)77. (31) 



z = 



Numerical estimates of this series based on the Pade- 
Borel resummation are given in the Table d for m = 1 
and compared with values found by other approaches. 

TABLE IV: The dynamical critical exponents of m-compo- 
nent model with extended impurities, obtained on the base of 
resumming the double e, Ed-expansion (Eq. i'27l ). 









m = 2 




m = 3 




m = 4 


Ed 




z± 


z \\ 


z± 


z \\ 


z± 


Z\ 





2.0145 




2.0150 




2.0151 




0.1 


2 


.021 


2.0147 


2.0150 


2.0150 


2.0151 


2.0151 


0.2 


2 


.046 


2.0150 


2.0150 


2.0150 


2.0151 


2.0151 


0.3 


2 


.067 


2.0151 


2.0511 


2.020 


2.0151 


2.0151 


0.4 


2 


.175 


2.094 


2.088 


2.050 


2.024 


1.984 


0.5 


2 


.209 


2.110 


2.129 


2.081 


2.064 


2.020 


0.6 


2 


.238 


2.122 


2.167 


2.101 


2.114 


2.054 


0.7 


2 


.270 


2.138 


2.210 


2.120 


2.158 


2.085 


0.8 


2 


.287 


2.146 


2.247 


2.140 


2.195 


2.113 


0.9 


2 


.315 


2.157 


2.282 


2.153 


2.244 


2.126 


1.0 


2 


.333 


2.171 


2.317 


2.169 


2.279 


2.141 


1.1 


2 


.354 


2.175 


2.346 


2.176 


2.314 


2.164 



V. EFFECTIVE CRITICAL BEHAVIOR 

The previous discussion concerns the asymptotic criti- 
cal dynamics of the model . Here we discuss the effec- 
tive critical behavior which is observed in approaching 
the critical point T c |2^|. This behavior is character- 
ized by effective critical exponents governing scaling laws 
when T c still is not reached. The calculation of effective 
critical exponents for models with extended impurities 
was not considered so far in dynamics as well as in the 
statics. However these exponents are mainly observed in 
numerical and real experiments, which may be outside 
the asymptotic region. 

The effective critical exponents are defined as logarith- 
mic derivatives of the appropriate quantities of interest 
with respect to the reduced temperature t = \T — T c \/T c 
|29|| . For instance, effective critical exponent for perpen- 



Substituting into (|31|l the e-expansion for rj 6lj, one gets: 



-0.363 0-^ + 0. 3 770\ g3 
(to + 8) 2 J 

I 

dicular correlation length is defined as: 

In the limit T — > T c the effective exponent coincides 
with the asymptotic one — v±. Within the field- 
theoretical RG approach effective exponents are calcu- 
lated using the values of couplings (solution of flow equa- 
tion) depending on the flow parameter £. For instance, 
effective exponents for correlation length and relaxation 
time perpendicular to the extended impurities direction 
are defined as: 

l/i/f(0 = 2-7*»(t*(*),«(0)- 
- l4> (u(£),v(£)) + ...; (34) 

zf(£) = 2 + ((u(£) 7 v(£)) + ... . (35) 

The flow parameter £ is related to the temperature dis- 
tance to the critical point via the inverse correlation 
length. Therefore the dependence of the effective ex- 
ponents on the flow parameter cor resp onds to a depen- 
dence on the reduced temperature [65[ . In l|3^|l and (f53f) 
the parts denoted by dots come from the change of the 
amplitude part of the perpendicular correlation length 
and critical relaxation time correspondingly. In the sub- 
sequent calculations we will neglect this part, since the 
contribution of the amplitude function to the crossover 
seems to be small [66) . 

Solving the flow equation for different initial conditions 
we can obtain different flows in the space of couplings. 
We choose these condition near the Gaussian FP (G), 
expecting that couplings are small in the background re- 
gion. In flow equations l|10l) we use the (3- function 120fl . 
(|21|l at fixed d resummed by the Chisholm-Borcl method. 
This allows us to investigate also the case of Ising systems 
(to = 1), for which the (3- function are degenerated on the 
one-loop level and the y^-expansion does not allow for 
numerical estimates |40j. 

The RG flows at to = 1 for the two most interesting 
cases Ed = and Ed = 1 are shown in the Fig. [21 (dashed 
lines for Ed — and solid lines for Ed — 1). The first case 
corresponds to the point defects, whereas the second one 
corresponds to the lines of defects. We consider the same 
initial condition for both cases with a small (flows 1 and 
1') and with a large (flows 2 and 2') value of the ratio 
v/u. The ratio v/u defines the degree of disorder, thus 



(TO + 2)e 2 / (to + 2)(-to 2 + 56to + 272) 

2 + 0.363 v ; „ + 0.09075 > M 



(to + 8) z \ (to + 8) 4 
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0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 2,0 2,2 2,4 2,6 2,8 




FIG. 2: Flows for the case m = 1. Dashed lines correspond to 
the case of uncorrelated point defects, Ed = 0, while solid lines 
correspond to system with extended lines of defects, Ed = 1- 
The initial values are the same for curves 1 and l' (small ratio 
A/w) and for 2 and 2' (large ratio w/u). 



we can observe difference in the behavior of systems with 
low and high dilution. 

Using the flows given in the Fig. we can obtain static 
and dynamic effective exponents. Below, we choose to 
display the exponents z/[" and z*f*. Due to Eq. I|19|) 
the parallel effective exponents i/fj and zjj have quali- 
tatively similar behavior. The deviations of effective ex- 
ponents from the mean- field values z c ^ — 2 and v± — 1/2 
corresponding to flows of Fig. [5] are given in the Fig. [3] 
and Fig. 0] respectively. First we compare the behavior 
of the effective exponents of the uniaxial m = 1 mag- 
nets at Ed = and Ed = 1. As it can be seen from the 
Figs. |3| and even a weak dilution by lines of defects 
leads to a faster crossover to the new universality class 
(curves 1 in Figs. |3and0J) and the effective exponents 
are not influenced by the pure FP P. Such a behavior 
qualitatively differs from that at dilution by point-like 
impurities, where effective exponents corresponding to 
the pure FP P are observed in a relatively wide region 
(curves 1' in both figures). Note that the dynamical crit- 
ical exponent at Ed — 1 approaches its asymptotics 
from above. Such a nonmonotonic behavior is typical in 
statics for system with point-like defects pi |4J|. How- 
ever we do not observe nonmonotonic behavior for 
It may be a feature of static effective critical exponents 
for system with extended impurities or consequence of 
our approximation. 

We also consider static and dynamic effective behavior 
for different values of order parameter dimension. Fig. [S] 
and Fig. present the z e ^ — 2 and i/[ ff — 1/2 correspond- 
ingly obtained for the same initial conditions and differ- 
ent values m — 1, 2, 3, 4. The nonmonotonic character 
of dependence of exponents on the logarithm of flow pa- 
rameter is observed only for z e ^ at m = 1. 



z-2 0.3 



-160 -140 -120 -100 



-60 -40 -20 



In / 



FIG. 3: Behavior of the dynamical effective exponent z"^ for 
m — 1 in the case of Ed — (dashed lines) and sj, = 1 (solid 
lines) . Curves correspond to flows of Fig|S] 



-1/2 




FIG. 4: Behavior of effective critical exponent of perpendic- 
ular correlation length v c ^ for m = 1 in the case of Ed = 
(dashed lines) and Ed — 1 (solid lines). Curves correspond to 
flows of FigH 



of magnetic systems with Ed-dimensional non-magnetic 
impurities of parallel orientation. The presence of such 
impurities leads to anisotropy in the systems; thus, two 
correlation lengths, parallel and perpendicular to the ex- 
tended defects exist, which diverge in the vicinity of crit- 
ical temperature with different exponents. The dynamic 



Z.-2 



m=1 



m=2 



m=3 



m=4 




In / 



VI. CONCLUSIONS 

The present study is dedicated to providing the reliable 
numerical estimates for the dynamical critical exponents 



FIG. 5: Behavior of dynamical effective exponent 2^ for 
system with extended impurities with Ed — 1 and different m. 
Initial conditions are the same in all cases (small ratio v/u). 
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m=2 




-100 -80 -60 -40 -20 

In I 

FIG. 6: Behavior of effective critical exponent of perpendicu- 
lar correlation lenght v±* for system with extended impurities 
with ed = 1 and different m. Initial conditions are the same 
in all cases (small ratio v/u). 

critical behavior is modified as well, so that two times of 
critical relaxation in directions parallel and transverse to 
extended impurities appear. The former results obtained 
for such systems are based on the double expansion in 
parameters e, Ed and thus are rather of a qualitative char- 
acter. 

An example of related systems is given by magnets 
with extended randomly oriented defects: these can be 
described by the model of Weinrib and Halperin [23| . 
There, the distribution of defects is characterized by the 
pair correlation function, falling off with distance for 
large separations according to a power law. The influence 
of such correlated defects on criticality was a subject of 
several theoretical studies |2^, an d predicted new uni- 
versal critical behavior found recently its confirmation in 
the MC simulations |2S| . 

We applied the resummation to the RG functions of 
the model, obtained in the minimal subtraction scheme 
EJ 113 ■ treating them directly for fixed d — 3 and 
fixed parameter £ d - The case Ed = describes point- 
like quenched disorder and reproduces well-known re- 
sults. For Ed > it was found, that the relation z± > z\\ 
holds for every Ed and m, which is connected to the fact, 
that the extended defects cut interacting paths of spins 



perpendicular to the extended-defect direction, so in the 
parallel direction the fluctuations are stronger. 

The data of Tables [I] IIVI give numerical values of the 
dynamical critical exponents for the model A dynamics 
of the TO-vector magnets with extended Ed-dimensional 
defects. Together with possible scenarios of the effective 
critical behavior discussed in the Section it should fa- 
cilitate experimental studies of the influence of extended 
defects on criticality. As we already noted, the point 
defects change the universality class at d — 3 only for 
the Ising magnets (m = 1). Values of the critical ex- 
ponent z of the pure and the diluted (by point-like im- 
purities) Ising model are given in Tables [n] IIIII respec- 
tively. Comparing the data in these tables one sees an 
increase of the exponent for the diluted model with re- 
spect to the pure one. This corresponds to a stronger 
divergency of the relaxation time (increase of the criti- 
cal slowing-down effects). Further, comparing Tables InTl 
and [I] one sees that extended impurities make this effect 
more pronounced, leading to further increase of both ex- 
ponents Zx,Zi\. However this increase (of the order of 8 
% if one compares z-exponents for point and line defects) 
is not dramatic to make MC simulations impossible due 
to enhanced critical slowing down. On the other hand, 
the same order of difference in critical exponents values 
holds also for the static exponents [2fll2ll |. Therefore, we 
expect that the predicted change in the critical behav- 
ior of three dimensional magnets with extended impuri- 
ties is within current experimental accuracy. Although it 
is more difficult to extract dynamical critical exponents 
from simulations to the same accuracy, it would be worth- 
while to look for them (i) in the cases of different defect 
dimensionality and (ii) in the non-asymptotic regime. 
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